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ABSTRACT 

We explore the effects of dark matter and dark energy on the dynamical scaling prop- 
erties of galaxy clusters. We investigate the cluster Faber- Jackson (FJ), Kormendy 
and Fundamental Plane (FP) relations between the mass, radius and velocity disper- 
sion of cluster size halos in cosmological A'^-body simulations. The simulations span a 
wide range of cosmological parameters, representing open, flat and closed Universes. 

Independently of the cosmology, we find that the simulated clusters are close to a 
perfect virial state and do indeed define a Fundamental Plane. The fitted parameters 
of the FJ, Kormendy and FP relationships do not show any significant dependence on 
r^TO and/or JIa- The one outstanding effect is the influence of Vim on the thickness of 
the Fundamental Plane. 

Following the time evolution of our models, we find slight changes of FJ and Ko- 
rmendy parameters in high Vim universe, along with a slight decrease of FP fitting 
parameters. We also see an initial increase of the FP thickness followed by a conver- 
gence to a nearly constant value. The epoch of convergence is later for higher values 
of Vim while the thickness remains constant in the low Vim A-models. We also find a 
continuous increase of the FP thickness in the Standard CDM (SCDM) cosmology. 
There is no evidence that these differences are due to the different power spectrum 
slope at cluster scales. 

From the point of view of the FP, there is little difference between clusters that 
quietly accreted their mass and those that underwent massive mergers. The principal 
effect of strong mergers is to change significantly the ratio of the half-mass radius Thai / 
to the harmonic mean radius Vh. 

Key words: Cosmology: theory - cosmological parameters - dark matter - large-scale 
structure of Universe - galaxies: clusters: general 



1 INTRODUCTION 

Recent observations of distant supernovae IjRiess et al ]|l998l : 
iPerlmutter et all 1 19991 ) suggest that we are living in a flat, 
accelerated Universe with a low matter density. This accel- 
erated expansion has established the possibility of a dark 
energy component which behaves like Einstein's cosmologi- 
cal constant A. A positive cosmological constant resolves the 
apparent conflict suggested by th e old age of globular clus - 
ter stars and the estimated value (jSpergel et al.|[2003l , |2007| ) 
appears sufficient to yield a flat geometry of our Universe. 

The role of A in the process of structure formation is 
not yet fully understood. Although its influence can be seen 
when looking at the global evolution of the Universe, its role 
in the dynamical evolution of cosmic structures is not clear. 
The most direct impact of A comes from its influence on the 
amplitude of the primordial perturbation power spectrum; 
there is also an influence from the change in the cosmic 
and dynamical time scales. The direct dynamical influence 



is probably minor: we do know that in the linear regime 
it accounts fo r a mere ~ l/70th of the influence of matter 
perturbation l|Lahav et al.lll99ll ). 

Most viable theories of cosmic structure formation in- 
volve hierarchical clustering. Small structures form first and 
they merge to give birth to bigger ones. The rate and history 
of this process is highly dependent on the amount of (dark) 
matter present in the Universe. In Universes with a low Vim, 
structure formation ceases at much early times than that in 
cosmologies with high density values. 

Within this hierarchical process, clusters of galaxies are 
the most massive and most recently formed structures in the 
Universe. Their collapse time is comparable to the age of the 
Universe. This makes them important probes for the study 
of cosmic structure formation and evolution. The hierarchi- 
cal clustering history from which galaxy clusters emerge in- 
volves a highly complex process of merging, accretion and 
virialization. In this paper we investigate in how far we can 
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get insight into this history on the basis of the internal prop- 
erties of the clusters. This involves characteristics like their 
mass and mass distribution, their size and their kinetic and 
gravitational potential energy. In particular, we are keen to 
learn whether these do show any possible trace of a cosmo- 
logical constant. 

One particular profound manifestation of the virial 
state of cosmic objects is via scaling relations that con- 
nect various structural properties. Scaling relations of col- 
lapsed and virialized objects relate two or three fundamental 
characteristics. The first involves a quantity measuring the 
amount of mass M, often expressed in terms of the amount 
of light L emitted by the object. The second quantity in- 
volves the size of the object, while the third one quantifies 
its dynamical state. For a virialized halo with mass M, size 
R and velocity dispersion ct„ = (u^)^''^, the implied scaling 
relation is 

log A'/ = 21og(7„ + logR + EM , (1) 

where ea/ is a constant that refiects the internal dynam- 
ics of the system. (eM is determined by issues such as the 
isotropy if the cluster velocity dispersion, its shape and any 
substructure) . 

Systems having similar values of this constant would be 
expected to form a two-parameter family of objects: obser- 
vationally this manifests itself as the "Fundamental Plane" . 
Objects lying on the same plane might be expected to have 
similar formation histories and, conversely, the nature of the 
Fundamental Plane is a clue to the underlying formation 
mechanism. 

The scaling relations are of great importance for a vari- 
ety of reasons. First of all, they inform us about the dynam- 
ical state of the objects and m ust be a profound reflec tion of 
the galaxy formation process ('Robertson et al]|2006l ). Also, 
they have turned out to be of substantial practical impor- 
tance. Because they relate an intrinsic distance independent 
quantity like velocity dispersion to a distance dependent one 
like Le, they can be used as cosmological distance indicators. 

1.1 Observed relationships 

1.1.1 Galaxies 

Since the mid 70s, we know that the observed properties of 
elliptical galaxies follow scaling re lations. The Faber-Jackson 
relation (|Faber fc JacksonI Il976l ) relates the luminosity L 
and the velocity dis persion a o f an elliptical galaxy. The 
TuUy-Fisher relation (|Tu11v fc Fi sher 1977) is the equivalent 
for spiral galaxies. A different, though related, scaling is that 
between the effective radius re and the luminosi ty L of the 
galax y. This is known as the Kormendy relation IjKormendvl 
1 19771 ). These two relations turned out to be manifestations 
of a deeper scaling relation between three fundamental char- 
acteristics, which became known as the Fundam ental Plane 
l|Diorgovski fc Davis|[l987l : iDressler et al]|l987l ) . 

The Fundamental Plane is generally expressed as a rela- 
tionship between three parameters, though there is no con- 
sensus as to which three should best be used, nor precisely 
how to deflne them. This makes detailed comparisons some- 
what diSicult. Some authors use the set (log i?, log cr, log 7) , 
I being the luminosity in some spectral band within some 
radius R, while others use the set (log _R, log ct, /i), fi being 



the mean surface brightness within that radius. Comparisons 
are further complicated by the fact that there appear to be 
manifest residual luminosity dependences in the fits, as re- 
ported in a recent study of the SDSS bv ,Nigoche-Netro et al.l 
j2009h . 

Care is needed when interpreting these observed re- 
lationships. Observed data generally refers to luminosity 
rather than mass, and the radius that is used generally refers 
to some fiducial radius such as the half-light radius or some 
radius based on profile fitting. Often, the half-light radius. 
Re as determined from a fit to a de Vaucouleurs profile is 
used. 

This situation has bee n improved s omewhat by the 
gravitational lensing study of lBolton et al.l |2007|). These au- 
thors presented a new formulation of the FP using lensing 
data to replace surface brightness with surface mass density. 
They also present an interesting alternative, which they refer 
to as the "Mass Plane" (MP), in which they find the depen- 
dence of log(_Re) on log(cre2) and surface mass density Ee2 
within a radius Re/2. Using surface mass density Ee2 within 
a radius Re/ 2 in place of surface brightness le removes one 
of the assumptions about the relationship between mass and 
light. 

1.1.2 Galaxy Clusters 

Much recent galaxy cluster work on the Fundamental Plane 
has focussed on the differences between the Fundamental 
Planes of the cluste rs as defined by their member galax- 
ies (see for example iD'Onofrio et al.l ( |200^ and references 
therein). 

Galaxy cluster s caling relations were discovered by 
ISchaeffer et all (|l993h who studied a sample of 16 galaxy 
clusters, concluding that these systems also populate a Fun- 
damental Plane. Adami et al. (1998) used the ESO Nearby 
Abell Cluster Survey (ENACS) to study the existence of 
a Fundamental Plane for rich galaxy clusters, finding that 
it is significantly di fferent from that for elliptical galaxies. 
iMarmo et al] (120041 ) using data from the WINGS cluster 
survey found that the difference is largely a simple shift in 
the relative positions of the planes. 

The largely unknown relationship between mass and 
light frustrates a direct comparison with the results of N- 
Body investigations. 

1.2 Numerical investigations 

Later, iLanzoni et al.l (|2004l ) addressed the question using 
N-Body simulations for high mass halos, which are thought 
to host clusters of galaxies. On the basis of 13 simulated 
massive dark matter halos in a ACDM cosmology they found 
that the dark matter halos follow the FJ, Kormendy and 
FP-like relations. 

In hierarchical scenarios of structure formation halos 
build up by subsequent merging of smaller halos into larger 
and larger halos. Some of these mergers involves sizeable 
clumps, most involves a more quiescent accretion of mat- 
ter and small clumps from the surroundings. This process 
leaves its mark on the phase-space structure of the halos. 
Indeed, these dark halo streams are a major source of atten- 
tion in present day studies of the formation of our Galaxy 
(jHelmi fc Whitelll999l : lHelmill2000l ). 
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Table 1. Cosmological parameters for the runs. The first column gives the identification of the runs, and the following columns give the 
present matter density parameter, the density parameter associated with the cosmological constant, fij. = 1 — Qm, — quantifies the 
curvature of the Universe, the age of the Universe in Gyr since the Big Bang, the mass per particle in units of lO^"/i~^M0, the mass 
cut of the groups given by HOP in units of lO^^/t^^MQ, the value of the (over)density needed to have virialized objects with respect to 
the background density, and similarly, but now with respect to the critical density. 



It remains an interesting question as to whether we can 
find evide nce for these merginR events in the Fundamen- 
tal Plane. iGonzalez-Garcfa fc van Albadal (|2003l ) look into 
the effects of major mergers on the Fundamental Plane and 
found that the Fundamental Plane does remain largely in- 
tact in the case of two merging ellipticals. However, what 
the effects will be of an incessant bombardment of a halo by 
material in its surroundings has not been studied in much 
detail. Given that this is a sensitive function of the cosmolog- 
ical scenario, we will study the influence on FP parameters 
and thickness in more detail. 

In this paper we address the specific question as to 
whether we can trace an infiuence of cosmic parameters 
in the scaling relations for simulated clusters, and in par- 
ticular the infiuence of the cosmic density parameter Qrn 
and the cosmological constant A. We use a set of dissipa- 
tionless A''-body simulations involving open, flat and closed 
Universes. All the simulations are variants of the cold dark 
matter (CDM) scenario, representing different cosmologies, 
concerning both different values for the mass density flm, 
for dark energy Qa and for the implied power spectrum of 
density perturbations and the related merging and accretion 
history of the clusters. 

The organization of this paper is as follows. In section 
[2] we describe the simulations and the deflnitions of the var- 
ious parameters we use. In section |3] we present a general 
description of the scaling relations which we investigate in 
this study before specifying the way in which we analyze 
them from the cluster-sized halos in our simulation. We in- 
vestigate the scaling relations of galaxy clusters in different 
cosmologies at z = in section |4] Section [5] addresses the 
evolution of the scaling relations as a function of redshift and 
cosmic time. We also investigate the dependence of merging 
and accretion on the scaling relations, which we discuss in 
section |6] The interpretation of our results on the Funda- 
mental Plane within the context of the virial theorem is 
discussed in section [T] Conclusions are presented in section 

m 



2 THE SIMULATIONS 

We perform thirteen N-body simulations that follows the 
dynamics of = 256^ particles in a periodic box of size 
L = 200/i~^Mpc. The initial conditions are generated with 
identical phases for Fourier components of the Gaussian ran- 
dom field. In this way, each cosmological model contains the 
same morphological structures. For all models we chose the 
same Hubble parameter, h = 0.7, and the same normaliza- 
tion of the power spectrum, ag — 0.8. The principal differ- 
ences between the simulations are the values of the mat- 
ter density and vacuum energy density parameters, flm and 
Qa- By combining these parameters, we get models describ- 
ing the three possible geometries of the Universe: open, flat 
and closed. The effect of having the same Hubble parameter 
and different cosmological constants translates into having 
different cosmic times. Table [1] lists the values of the cosmo- 
logical parameters and the cosmic times at which the data 
is analysed. 

The initial conditions are evolved up to the present time 
{z = 0) using the mas sive parallel tree N-body code GAD- 
GET2 l|Springel|[2005l) . The Plummer-equivalent softening 
was set at tpi = 15ft~^kpc in physical units from z = 2 to 
z — 0, while it was taken to be fixed in comoving units at 
higher redshifts. For each cosmological model we wrote the 
output of 100 snapshots, from aexp = 0.2 {z = 4) to the 
present time, aexp = 1 [z — 0), equally spaced in log(a). 

2.1 Halo identification 

We use the HOP algorithm (|Eisenstein fc Hut|[l99^ ') to ex- 
tract the groups present in the simulations. HOP associates 
a density to every particle. In a first step, a group is de- 
fined as a collection of particles linked to a local density 
maximum. To make a distinction between a high density 
region and its surroundings, HOP uses a regrouping pro- 
cedure. This procedure identifies a group as an individual 
object on the basis of a specific density value. For this crit- 
ical value we chose the virial density value Ac based on the 
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Figure 1. Cluster evolution: Qm influence. Evolution as a function of redshift of a single dark matter halo in diflferent cosmological 
models: ACDM02, ACDMF2, ACDMC2 and SCDM. Shown is the dark matter particle distribution in a box of comoving size 5h~^ 
Mpc, at 6 different redshifts: z=2.98, z=2.28, z=1.53, z=1.01, z=0.50 and z=0.00. The circles correspond to halos identified by HOP, 
with the size of the circle being proportional to their virial radius. 
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Table 2. Cosmic times in Gyr and the corresponding redshift for 
a set of four reference cosmological models. 



spherical collapse model. In order to have the proper Ac we 
numerically compute its value for each of the cosmologies. 
See Table [T] for the values of the virial density for each cos- 
mology at 2; = 0. For the latter we list two values: the virial 
overdensity A„ir,6 with respect to the background density 
Pb of the corresponding cosmology, and the related virial 
overdensity Avir.c with respect to the critical density. 

Note that we only consider groups containing more than 
100 particles. Because the particle mass depends on the cos- 
mological scenario, this implies a different mass cut for the 
halos in each of our simulations. As a result, SCDM does 
not have groups with masses lower than 10^'^h~^MQ. We 
have to keep in mind this artificial constraint when consid- 
ering collapse and virialization in hierarchical scenarios at 
high redshifts, and also when making fits to the relation- 
ships among the various cluster parameters. In cases where 
structure growth is still continuing vigorously at the current 
epoch, the collapsed halos at high redshifts will have been 
small; our simulations would not be able to resolve these. 



2.2 Halos and Cosmology: an example 

Figure [1] follows the evolution of one particular cluster 
halo in four different cosmologies. These are ACDM02, 
ACDMF2, ACDMC2 and SCDM. By using the same Fourier 
phases to set up the initial conditions in each of the cosmolo- 
gies we get a sample of corresponding halos. In each of the 
cosmologies the evolution of the cluster halo is shown at six 
different redshifts, from 2: r; 3 onward to the present epoch 
z = 0. The panels show the mass distribution in and around 
the cluster, and its progenitors, in a box of comoving size 
5h~^ Mpc. Circles enclose halos identified by HOP, with the 
circle radius proportional to the virial radius of the group 
(i.e. the distance from the center of mass to the outermost 
particle of the group). Projection effects may occasionally 
cause circles to appear within circles. 

In all four cosmologies, the buildup of the halo clearly 
involves the merging of several smaller mass clumps, some 
of which are identified as genuine proto-halos by means of 
circles. Fig.[T]shows that the sequence ACDM02, ACDMF2, 
ACDMC2 and SCDM corresponds to a sequence in which 
the formation of the halo shifts to later and later epochs. At 
all depicted redshifts, and in particular at higher redshifts, 
the clusters in the ACDM02 cosmology have considerably 
more pronounced and developed mass concentration. 



2.3 Halo properties 

In our study, we limit ourselves to cluster-like halos. A 
galaxy cluster is defined as a dark matter halo with a mass 
M> lO^*/i~^M0. We measure three quantities for each clus- 
ter and test their scaling relations. 

Scaling relations of collapsed and virialized objects re- 
late two or three fundamental characteristics of those ob- 
jects. The first involves a quantity measuring the amount of 
mass, often in terms of the amount of light emitted by the 
object. The second quantity involves the size of the object, 
while the third one quantifies its dynamical state. 

• Mass: defined as the number of particles multiplied by 
the mass per particle present in each group: 



M 



(2) 



where Upart is the number of particles in the halo and nipart 
is the mass of each particle (see column in Table 

The mass of the particle is different for each cosmology. 

• Surface mass density: Alternatively, following obser- 
vational practice, we use the magnitude-scale surface mass 
density /i for our Fundamental Plane evaluations. 



^ = -2.5 log A'/ -I- 51ogr 



(3) 



where M and r are the mass and the radius of the halo. 
Combining this with a mass to light ratio it becomes a sur- 
face brightness, one of the observables of the Fundamental 
Plane. 



Velocity dispersion: computed as 
2K 



f^part TTT'part 

where K is the kinetic energy of the halo. 



(4) 



As a measure for the size of the halos, we have explored two 
options: the half-mass radius and the mean harmonic radius. 

• Half-mass radius: rhaif is the radius that encloses 
half of the mass of the clump. This radius is closest in 
definition to the half-light radius used in observational 
studies. 

• Mean harmonic radius: rh is defined as the inverse 
of the mean distance between all pairs of particles in the 
halo: 



1 



-y— 

i<3 



N = 



f^parti^f^npart 



(5) 



where r^j is the separation vector between the ith and the 
jth particle. The great virtue of this radius is that it is a good 
measure of the effective radius of the gravitational potential 
of the clump, certainly important when assessing the virial 
status of the clump. Also, it has the practical advantage of 
being independent of the definition of the cluster center. To 
some extent, it is also an indicator of the internal structure 
of the halo because it put extra weight to close pairs of 
particles. 

Most of the results presented in this paper refer to the 
mean harmonic radius of the halos: this seems rather natural 
given that we are discussing the virial theorem (see Table 
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3 SCALING RELATIONS 
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Figure 2. Comparison between the mean harmonic radius and 
the half-mass radius of the cluster-size halos in the SCDM (top) 
and ACDMF2 (bottom) scenarios. The colours depict different 
mass ranges, each colour representing a 20 percentile mass quan- 
tile. 



[3]). We have also compared the results obtained using the 
half-mass radii of the halos. 

In Fig. [2] we plot the mean harmonic radius versus the 
half-mass radius of the cluster-size halos in the SCDM and 
in the ACDMF2 models. We see that the relationship is not 
very tight at larger masses, and that the differences between 
the two radii are particularly prominent in the SCDM cos- 
mology. We shall discuss this further in section [S] Not sur- 
prisingly, the fitted Fundamental Plane parameters depends 
strongly on which radius is used. Equally surprising, the Ko- 
rmendy relation slope does not seem to be particularly sen- 
sitive to the choice of or r^ai / (the slopes are statistically 
not different). This is summarised in Table [G] 



For the samples of cluster-sized halos in our simulations we 
will be assessing three specific scaling relations: the Faber- 
Jackson relation, the Kormendy relation and the Fundamen- 
tal Plane. 

From observations of elliptical galaxies we have 
learned that there are tight scaling relations between a 
few of their fundamental structural properties (see e.g. 
iBinnev &: Merrifieldl (|l998 l)). These properties are the to- 
tal luminosity L of a galaxy - or its surface brightness - 
its characteristic size Re and its velocity dispersion o-„. 



3.1 Faber- Jackson and Kormendy Relations 

The first scaling relat ion is the Faber- Jackson relation 
(|Faber fc JacksonI \l97d ) between the luminosity L of the 
galaxy and its velocity dispersion (t„. 



L oc (t: 



/3 



(6) 



where the index /3 4. A similar relation , known as the 
TuUy-Fisher relation (ITuIIv fc Fisheij Il977^ ■ holds for HI 
disks of spiral galaxies. According to this relation, the galax- 
ies' rotation velocity is tightly correlated with the absolute 

magnitude of the galaxy. 

Another relation was established bv lKormendvl l|l97'i1 ). 
He found that there is a strong, not entirely unexpected, 
correlation between the luminosity L and effective radius 
Re of the elliptical galaxies: 

Lo^Re , (7) 
where the index a ~ 1.5. 

3.2 Galaxy Fundamental Plane 

Both the FJ and Kormendy relations relate two structural 
characteristics and should be seen as projections of a more 
fundamental and tight relation between all three struc- 
tural quantities: the Fundamental Plane (FP). The Funda- 
mental Plane of elliptical gala xies was first formul ated by 
iDiorgovski fc Davis! l|l987l ) and lDressler et al.l (|l987D . When 
we take the three-dimensional space defined by the effective 
radius Re of the galaxy, its surface brightness le (with to- 
tal luminosity L oc leRe) and velocity dispersion ct„, we find 
that they do not fill space homogeneously but instead define 
a thin plane. 

In logarithmic quantities, this plane may be parameter- 
ized as 

log -Re = 7log-^e + Sloga^ + Cfp (8) 

For example, |j0rgensen et al] l| 19961 ) found that a reasonable 
fit to the Fundamental Plane is given by 

logi?e = -0.82 log /e + 1.24 log a„ + Cfp (9) 

While nearly all galaxies, ranging from giant ellipti- 
cals to compa ct dwarf ellipticals, ap pear to lie on the FP 
(also see e.g. j0rgensen et al. 1995; Bcrn ardi et al. I l2003l : 
ICappellari et al. 2006: Bolt on et al. 2007) it is interesting t o 
note that diffuse dwarf ellipticals do not (lKormendvlll987l ): 
they seem to be fundamentally different objects. 

The observed Fundamental Plane not only provides in- 
formation on the dynamical state of the object but also on 
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1.69 


7.93 


0.01 


6.97 


ACDM03 


0.3 


0.5 


0.36 


0.05 


0.34 


0.03 


0.41 


1.03 


1.86 


3.88 


0.03 


11.43 


ACDMF2 


0.3 


0.7 


0.36 


0.05 


0.33 


0.03 


0.41 


1.11 


1.88 


4.01 


0.02 


11.23 


ACDMCl 


0.3 


0.9 


0.34 


0.05 


0.34 


0.03 


0.42 


1.19 


1.92 


4.33 


0.03 


11.09 


ACDMF3 


0.5 


0.5 


0.35 


0.05 


0.34 


0.03 


0.38 


0.56 


1.81 


2.10 


0.03 


11.75 


ACDMC2 


0.5 


0.7 


0.35 


0.05 


0.34 


0.03 


0.38 


0.57 


1.82 


2.14 


0.03 


11.66 


ACDMC3 


0.5 


0.9 


0.34 


0.05 


0.34 


0.03 


0.39 


0.61 


1.83 


2.31 


0.03 


12.24 



Table 3. Parameters of the scaling relations derived for the galaxy clusters in each of the simulated cosmological simulations, is the 
mean harmonic radius of the cluster, a is the scaling parameter for the Kormendy relation, b is the scaling parameter for the Faber- 
Jackson relation and c and d are the scaling parameters for the Fundamental Plane, a is the standard error in each of the scaling relation 
parameters, S is the corresponding standard error/significance of the fit. 



the evolution of its stellar content and, by implication, about 
its formation. For a virialized object with effective radius Re 
and mass-to-light ratio M/L the FP relation will have the 
form 

logiJe = - log/, + 2 log (J, - log(A//L) + Cs , (10) 

in which le = L/A-kR^ is the mean surface brightness and 
Ca a constant dependent on the structure of the object. 

The observed parameter values for elliptical galaxy Fun- 
damental Plane (see Eqn. [UJ are different from what might 
be expected for a plane that results simply from virializa- 
tion and constant mass-to-light ratio. One explanation for 
this difference is that galaxies may be structurally equiva- 
lent while having a mass- dependent M/L ratio. That would 
imply a formation process involving a tight fine tuning of 
M/L. Ne vertheless, pursu i ng thi s view, the parameters in- 
ferred by |j0rgensen et al.l (|l996l ) (Eqn. |9]) would imply a 
mass-to-light ratio dependence on mass: 

(M/L) (X M°-^'' , (11) 

using M (x a^Re and L oc leR^ (see e.g. iFabeij Il987l ). 
Recent semi-analytical modelling of galaxy formation sug- 
gest a more complex relation between the mass-to- light ratio 
and luminosity, involving a minimum M/L for galaxies with 
M « lO" - 10^2/j-iMg^ the absence of any mass-to-light 
dependency, the discrepancy between the planes would have 
to be due to variations in the structure parameters of the 
galaxies. 

There is an intrinsic scatter of the FP that has been 
found for elliptical galaxies: this has not been completely 
explained and may be a manifestation of the formation pro- 
cess. 

A slightly dif ferent approach is us ed in the gravitational 
lensing study of iBolton et al.l (jiooj) • These authors pre- 
sented a new formulation of the FP using lensing data to 
replace surface brightness with surface mass density, arriv- 
ing at the relationship of the form 



log Re = 7log/e + Slogae2 + Cfp, (12) 

where (Te2 is the velocity dispersion within half of the effec- 
tive radius Re, and 

7 = -0.78 ±0.13, 5 = 1.50 ±0.32, C/p = 3.9 ± 1.7. (13) 

Furthermore they suggest that the scatter about the Funda- 
mental Plane, derived from their data, correlates with their 
derived mass-to-light ratio for the galaxies in their sample. 
The evidence is not strong though it is suggestive. 

They also present an interesting alternative, which they 
refer to as the "Mass Plane" (MP), in which they find the 
dependence of log{Re) on log(CTe2) and surface mass density 
Ee2 within a radius Re/2: 

log Re = 7mlogEe2 ± <5m log Cre2 ± Cfp,,n (14) 

with 

7„ = -1.16 ± 0.09, 5„ = 1.77 ± 0.14, C/p = 7.8 ± 1.0. (15) 

Using surface mass density Ee2 within a radius Re/ 2 in 
place of surface brightness Le removes one of the assump- 
tions about the relationship between mass and light. 

3.3 Cluster Fundamental Plane 

If clusters were fully virialized objects with the same internal 
dynamics, they would necessarily lie on a universal Funda- 
mental Plane in t he Mass-velocit y -radiu s space. This was 
first addressed by ISchaeffer etabl (Il993l ). who, using sam- 
ple of 29 Abell clusters, discovered a FP relation in light- 
velocity-radius space : L oc Re'^^ al'^^ . This is equivalent to 
the relationship 

logi?e = -0.90 log ± 1.15 log cr„ ± C/p,, (16) 

in which le is a measure of the mean surface brightness of 
the cluster. The corresponding FJ relation is L oc and 
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the Korme ndy relation is L oc R l '^"^. Similar numbers were 

inferred bv lLanzoni et al.l l|2004 ). L oc i?°-^°(Tif^ 

In a project designed to test this further, lAd ami et alJ 

found a FP relation for a sample of ENACS Clusters, 
though their fitted parameters were markedly different: L oc 
^i.i9±o.i4^o.9i±o.i6^ -pjjjg jg equivalent to the relationship 



log Re 



(1.23±0.20) log Je + (1.12±0.11) logcr„ +C/P , (f7) 



in which is the mean surface brightness of the cluster. Note 
that there are considerable systematic uncertainties in these 
values which are not reflected in the quoted error bars: these 
arise out of the profile fitting to the cluster. The above fit 
to the data is based on fitting a King profile, (this gave the 
best fit to the data). 

In studies of simulated dark matter dominated galaxy 
clusters, we can study scaling relations that are similar to 
those inferred from observable quantities. To infer these rela- 
tions we base ourselves on the mass M of the object. If the 
selected objects have the same average density, we would 
expect an equivalent Kormendy relation given by 



M (xRt 



(18) 



Any difference in slope should be ascribed to a dependence of 
mean density {p{Re)) on the size Re of the object. The equiv- 
alent Fundamental Plane relation will be that of Eqn. [Tl 
while the Faber- Jackson relation would then be 



M oc at 



(19) 



Note that this is based on the assumption of constant mean 
density p of the selected objects, in line with HOP overden- 
sity criterion (see sec. 12.111 . 

iLanzoni et al.l l|2004) analyzed the N-Body cluster scal- 
ing relations on the basis of a sample of 13 massive dark 
matter halos identified in a high resolution ACDM A'^-body 
simulations. They were able to confirm the existence of FP 
relations for the simulation dark matter clusters and also 
found that these have a slope that was significantly differ- 
ent from the galaxy FP slope, 

logRe = (0.44 ± 0.02) 4- (1.92 ± 0.12) log cr„ + C/p , (20) 

with /i the surface mass density (Eqn. [5}. The difference 
in FP parameters between the dark matter halos and those 
inferred for the observed cluster sample (see above, Eqn. ll7p 
formed a key aspect of their study. They suggest a mass 
dependent cluster M/L ratio 



(M/L) oc M° 



(21) 



would be able to explain the observed cluster Fundamental 
Plane. Interestingly, this is markedly different from that in- 
ferred for early type galaxies. Of course there is no obvious 
reason why the Fundamental Plane for galaxies should have 
any bearing on the Fundamental Plane for clusters. Indeed, 
as we shall see for the ENACS sample, its FP parameters 
values seem to be irreconcilable with the virial theorem. 



3.4 Determination of Scaling Relations 

For the sample of N cluster-sized halos in each simulation we 
study the scaling relations between their size r, mass M - or 
equivalent surface mass density p - and velocity dispersion 
a (note that N is in general different for each cosmology). 



Given the inferred mass M (Eqn. [2)| , velocity dispersion (7„ 
(Eqn. |4]) and the mean harmonic radius rn (Eqn. [5]) of the 
cluster halos, we find the scaling relation parameters by lin- 
ear fitting of the relations. 

Sample selection effec ts play a complex r ole in the anal- 
ysis of real data samples (|La Barbera et al] (|2003l )'). Fortu- 
nately the issue is far simpler when analysing clusters found 
in N-Body models where the only selection criterion is a 
mass cut-off imposed by the cluster finding algorithm. We 
deal with that simply by making the Mass of the cluster the 
independent variable in all fits where relevant: this elimi- 
nates biases introduced through this object selection. 



3.4-1 Kormendy Relation 

For the Kormendy relation we fit 

logr = a logA/ -I- C„. (22) 

to the N data points (log r^, log KU) of the halo sample. The 
significance Sk is computed from the A'^ residuals: 



Sk 



a\oghh~Ca)2 



(23) 



3.4-2 Faber- Jackson Relation 

Along the same line, the Faber- Jackson relation is deter- 
mined on the basis of the fit 



log a„ = 6 log M + Cb, 

whose significance iSfjIs calculated as follows: 



Sfj = 



JV 



\ (A- 



1) ^ 



(logCTt,,i - 6 log Mi - Cb)2 



(24) 



(25) 



3.4-3 Fundamental Plane 

Instead of fitting the Fundamental Plane in the form of 
Eqn. [1] we do it in the way suggested by the observational 
work, i.e., using the surface mass density /i and velocity dis- 
persion CTi, as free parameters from which we determine a 
model for the radius, 



log r = c ^ 4- d log cr„ -f- Cfp . 



(26) 



In this, /i is the magnitude-scale surface mass density 
(Eqn. [31). Although there are errors in determining both ct„ 
and /i, they are very small when compared with the disper- 
sion about the Fundamental Plane. By fitting the parame- 
ters c and d this way we solve problems regarding biases in 
the mass (luminosity) selection. 

The significance Sfp of the Fundamental Plane fits de- 
rived from the sample of N simulation cluster halos is com- 
puted according to: 



(iV-2) 



(log ri - cpi - log cr„ - C/p)2 . (27) 
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Figure 3. Kormendy Relation. Each panel plots the relation between mean harmonic radius rj, and mass M of the cluster-sized dark 
halos in the simulations corresponding to one particular cosmology. Going from top left to bottom right these are: SCDM, ACDM02, 
ACDMF2 and ACDMC2. In each of the panels we have superimposed the fitted Kormendy relation for the corresponding model and for 
ACDMF2 as comparison. 



The thickness w fp of the Fundamental Plane is estimated on 
the basis of the perpendicular distances of the cluster halos 
to the fitted plane: 



Wfp 



E3 

N 



(28) 



where TV is the number of cluster halos in the sample and 
D± is the perpendicular distance of a point to a plane 



Da 



cfi + d log ay + Cfp - log Th 

(c2 +d2 + 1)1/2 



(29) 



4 SCALING RELATIONS IN DIFFERENT 
COSMOLOGIES: z=0 

We first investigate the scaling relations of the cluster dark 
matter halos in our cosmological models at the current 
epoch, z — 0, and look for possible systematic differences 



between the parameter values and FP thickness as a func- 
tion of the cosmology. The parameters of the resulting linear 
fits, to be discussed in the following subsections, are listed 
in Table |3 



4.1 Kormendy Relation 

Fig.Oshows the relation between the mean harmonic radius 
Vh of each cluster halo and their mass M. Each of the four 
panels depicts the relation for the halos in one particular 
simulated cosmology. The top left panel shows the SCDM 
cosmology, the top right one the ACDM02 model, the bot- 
tom left one the ACDMF2 model and the bottom right one 
the ACDMC2 model. 

In each cosmology there is a strong and systematic 
almost linear relation between logM and logr^: the Kor- 
mendy relation appears to be a good description for all sit- 
uations. A visual comparison between SCDM relation (top 
left panel), the ACDM02 relation (top right panel) and the 
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Figure 4. Inferred parameter a for the Kormendy relation 
(Eqn. I22II as a function of Cm (top panel), Ha (central panel) 
and Om + f^A (bottom panel). The bars represent the la uncer- 
tainty range around the estimated parameter. 
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Figure 5. Inferred scaling parameter b for the FJ relation as a 
function of three different parameters: Vim (top panel), £7j\ (cen- 
tral panel) and Q,m + t^A (bottom panel). The bars represent the 
Icr uncertainty range around the estimated parameter. 



ACDMF2 relation (bottom left panel) shows that clusters of 
comparable mass have a larger size in the low Q.m cosmology 
than in the ones with a higher density value. In other words, 
clusters are more compact in the SCDM cosmology. Not un- 
expectedly we find objects of a higher density in higher Q,m 
models. 

When fitting the plotted point distributions, we infer 
the parameter values listed in Table |3] In each of the pan- 
els in Fig. [3] we plotted the linear fits for all of the four 
depicted cosmologies. We find similar slopes for all cosmolo- 
gies, in the order of a 0.36 — 0.38. This seems to imply 
that the mean density (p(r'h)) oc Af~"'^; more massive halos 
have a slightly lower average density (see also iLanzoni et al.l 
l|2004 )). To investigate the dependence of the Kormendy pa- 
rameter a on the cosmology in Fig. |4] we have plotted the 
slope o as a function of the average mass density parameter 
Q.m (top panel), as a function of the cosmological constant 
f^A (central panel) and as a function of the cosmic curvature. 



in terms of Q,totai ~ + (lower panel) . There is no evi- 
dence for any systematic trends of the Kormendy parameter 
as a function of cosmology. No evidence for an influence of 
either cosmic density Sim and JIa on the internal structure 
of the halos could be detected. 



4.2 Faber-Jackson Relation 

Fig. [6] shows the Faber-Jackson relation: the relation be- 
tween the mass M and the velocity dispersion (t„ of the 
cluster halos. Like in Fig. |3l each of the four panels corre- 
sponds to one particular simulated cosmology: SCDM (top 
left panel), ACDM02 (top right panel), ACDMF2 (bottom 
left panel) and ACDMC2 (bottom right panel). 

For comparison, in each of the panels we show the line 
of the ACDMF2 model corresponding to the linear fit of this 
relation in each of the depicted cosmologies. The M — re- 
lation is clearly well fitted by the Faber-Jackson like relation. 
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Figure 6. Faber-Jackson relation. Eacli panel plots the relation between the velocity dispersion ay and the mass M of the cluster-sized 
dark halos in the simulations corresponding to one particular cosmology. Going from top left to bottom right these are: SCDM, ACDM02, 
ACDMF2 and ACDMC2. In each of the panels we have superimposed the fitted Faber-Jackson relation for the corresponding model and 
for ACDMF2 for comparison. 



It is considerably tighter than the equivalent Kormendy re- 
lation. 

It is also interesting is to note that, as with the Kor- 
mendy relation, we do not find any significant dependence 
of the FJ relation on the underlying cosmology: the slope b 
in all cases is in the order of fe ~ 0.35 (see Table [Sj. We also 
did not find any dependence on Qa or iltotai (see Fig. 

Although the difference between the inferred value of 
b 0.35 in most cosmologies and the value of 6 = 0.33 
expected for virialized perfectly homologous systems (see 
Eqn. ll9|l is not really significant, the consistent and system- 
atic value b > 0.33 might be suggestive for a we akly homolo- 
gous p opulation along the lines described in e.g. lBertin et al.l 
(2002h . 

4.3 Fundamental Plane 

The Kormendy relation and the Faber-Jackson relation are 
two dimensional projections of an intrinsically three dimen- 



sional relation between mass M, size r and velocity disper- 
sion (T„ of the halos. By implication, the spread of the Fun- 
damental Plane relation should be less than that of each of 
the two previous relations. 

The Fundamental Plane obtained for the same cosmolo- 
gies as shown in Fig. [3] and |6] (SCDM, ACDM02, ACDMF2 
and ACDMC2) is illustrated in Fig. [7] In each of the frames 
we have plotted the harmonic radius of the halos against 
the quantity Y = cn + d log + Cf p on a log-log plot. The 
parameters c and d in the latter quantity, Y, combining the 
surface mass density ^ and the velocity dispersion ct„ of each 
halo, are the best fit FP parameters for the corresponding 
cosmology (see Table [S}. 

The galaxy clusters in each cosmology do indeed seem 
to populate a tightly defined plane. The point clouds in each 
of the frames confirm our expectation that they should have 
a much lower scatter around the plane than in the case of 
the Kormendy and Faber-Jackson relation (see Table [3]). 

From Table |3] we find a surprising level of consistency 
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Figure 7. Fundamental Plane. Each panel plots the relation between harmonic radius r^j and the quantity Y = cfi + d log cr. Combining 
the surface mass density fi and velocity dispersion ctv , the scaling parameters c and d are the ones inferred from the FP fitting procedure. 
Top left panel: the relation between harmonic halo radius and Y for the cluster halo sample in the SCDM simulation. Top right panel: 
for the halos in the ACDM02 simulation. Bottom left panel: for the halos in the ACDMF2 simulation. Bottom right panel: for the halos 
in the ACDMC2. The superimposed lines in each panel represent the relation for the fitted Fundamental Plane for the corresponding 
cosmology. Note that, by definition, each of these fitted lines should have slope unity. 



between the Fundamental Planes in each of the cosmolo- 
gies. We find that the inferred parameters are close to the 
one theoreticaUy expected for perfectly homologous virial- 
ized clusters halos. The inferred scaling parameter c for the 
surface density hovers around 0.38 — 0.42, close to the the- 
oretical value c ~ 0.4 (M oc Thcrt)- The difference is some- 
what larger for the parameter d, implying that the velocity 
dispersion scaling has a difference of ~ 0.15 — 0.25 from the 
theoretical value of 2. 

As can be seen in both Table [3] and Fig. [71 there is 
hardly any variation between the FP relations in the differ- 
ent cosmologies: they almost all coincide. This is certainly 
true concerning the FP parameters c and d. The two top 
panels of Fig. [5] do confirm the impression that there is no 
systematic difference as a function of fim and/or JIa- This in 
itself is a strong argument against differences in the scaling 



relations parameters being due to a p artial or incomplet e 
level of virialization, as was claimed bv iAdami et al.l (|l998h . 

One possible difference between the Fundamental Plane 
in different cosmologies may concern its thickness uijp. In- 
spection of Fig. [7] does suggest a marginally lower thickness 
of the FP for Universes with a low Q.m ~ 0.1. There is no 
detectable effect at all with respect to the cosmological con- 
stant Q,A- We might understand a dependence on ^m, or 
cosmological constant Qa, in terms of the ongoing evolution 
of the cluster population. In low Q,m Universes - and in high 
Q.A universes all clusters formed at high redshift and have 
since had ample time to reach full virialization and hence 
tighten the corresponding Fundamental Plane. In high Q.^ 
Universes, clusters would still undergo a substantial levels of 
merging and accretion, both of which may affect the virial 
state of the cluster. Our computer experiments do not seem 
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Figure 8. Fundamental Plane relation of dark halos in the SCDM (left) and ACDM02 (right) cosmology, using the half-mass radius 
'''half (top row) and the harmonic radius (bottom row). Plotted are mass M versus the FP quantity Y = + dlogcr, with c and d 
the scaling parameters inferred from the (linear) fitting procedure. The lines represents the best fit FP relations. 
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0.38 
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0.36 


0.08 
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0.1 
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Harmonic 


0.38 


0.05 


0.38 


1.66 


0.02 








Half-mass 


0.35 


0.07 


0.31 


1.66 


0.03 


ACDMF2 


0.3 


0.7 


Harmonic 


0.36 


0.05 


0.41 


1.88 


0.02 








Half-mass 


0.35 


0.08 


0.30 


1.63 


0.03 


ACDMC3 


0.5 


0.9 


Harmonic 


0.35 


0.05 


0.38 


1.82 


0.03 



Table 4. Scaling relation parameters and radius definition: inferred Kormendy relation parameter a and Fundamental Plane parameters 
c and d, based on the use of half mass radius Tf^^l f ^^nd harmonic radius r^. For four different cosmologies - SCDM, ACDM02, ACDMF2 
and ACDMC3 - the scaling parameters and the corresponding goodness-of-fit S (see Eon. I23l and Eon. I27fl are given for rhaif (top row) 
and rjj (bottom row). 
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Figure 9. Top panel: Fundamental Plane parameters c (left axis, 
solid circles) and d (righthand axis, solid triangles) as a function 
of f2m- Second panel: Fundamental Plane parameters c (left axis, 
solid circles) and d (righthand axis, solid triangles) as a function 
of Qa- Third panel: thickness wjp of the Fundamental Plane, i.e. 
rms scatter of the FP relation as a function of Om- Bottom panel: 
thickness Wfp of the Fundamental Plane, ie. rms scatter of the 
FP relation as a function of Ha- 



to find any strong and significant dependence on overall cos- 
mology. 

We investigate the relationsliip between the FP thick- 
ness and the dynamical state of the cluster in more detail in 
section [G] 

Finally, we can try to relate the Fundamental Plane (/i, 
T'h, of our simulated cluster samples to the observation- 
ally measured (L, Re, (t„) plane, e.g. L oc R^'^^a^'^^ found 
for the EN ACS survey.. We can ask whether the difference 
can be ascribed solely to a mass dependent mass-to-light 
ratio M/L. 



4.4 Scaling Relations for alternative Radius 
Definition 

Apart from the mean harmonic radius that we have used 
as a measure of halo size in the previous sections, we have 
also assessed the viability of the scaling relations in case of 
alternative size definitions. In Table we list the resulting 
parameters for the Kormendy relation and the Fundamental 
Plane in the case of using the half-mass radius rhaif- 

The parameters for the Kormendy relation hardly differ 
from the ones inferred on the basis of the mean harmonic 
radius. However, the inferred Fundamental Plane plane pa- 
rameters do difi'er significantly from the ones inferred above 
on the basis of the mean harmonic radius. The change in 
scaling parameter values may be ascribed to the use of quan- 
tities that probe different aspect of the structure and dynam- 
ics of the halos. In an extreme situation, this might have 
disrupted the scaling relations. Our finding shows that the 
Kormendy relation still holds, while the FP relation still 
holds but in a slightly different guise. It may be an indica- 
tion for our contention that halos do not form a perfectly 
a perfectly homologous population. Size measures sensitive 
to different aspects of the halos' internal mass distribution 
may then result in somewhat different scaling properties. In 
this respec t, we agree with the conclusions of lAdami et al.l 
(1998) and lLanzoni et all (|2004l ). 

See section [6] and [7] for a discussion of the relationship 
between the radii rhaif and r^, where we show that it is a 
consequence of the cluster building process. 



5 EVOLUTION OF SCALING RELATIONS 

In the previous sections we have extensively studied the scal- 
ing relations at the current cosmic epoch z = 0. We have 
also noted that there are differences between the scaling re- 
lation parameters that we find in our simulations and those 
for perfect virialized and homologous systems. This makes 
it interesting to trace the evolution of the different scaling 
relations. 

In this section we investigate the evolution of the scaling 
relations as a function of redshift and as a function of cos- 
mic look-back time. While observers usually think in terms 
of redshift, it is important to appreciate that a given red- 
shift corresponds to an entirely different dynamical epoch 
in different cosmologies. Given the same Hubble parameter, 
the age of the Universe is a sensitive function of the cosmic 
density parameter Qm and even more so of the cosmological 
constant. As for the latter, we have to realize that the change 
in cosmic time as a function of the cosmological constant is 
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Figure 10. Evolution of the fitted scaling relation parameters as a function of redshift (left column) and as a function of cosmic look-back 
time (right column). Top: Kormendy parameter a. Center: Faber-Jackson parameter b. Bottom: FP parameters c and d. 



16 Araya-Melo et al. 



0.015 




ACDM02 
ACDMF2 
ACDMC2 



0.000 , , , , , , , I 

0.2 0.4 0.6 0.8 1.0 

^exp 



Figure 11. Evolution of the thickness of the Fundamental Plane for four different cosmologies. Note the almost consistently tighter FP 
for the low Qrn Universe and the modest increase of FP thickness in the other cosmologies. The SCDM, ACDMF2 and ACDMC2 models 
are very similar in their behaviour. The more erratic behaviour of the ACDM02 may in part be due to the smaller sample size 



the most important influence of A. To give an appreciation 
of the differences in cosmic time for a given redshift in the 
different cosmologies, we refer to Tabled 

We have probed the scaling relations over a range of 
redshifts from z — 4 to z = and over a range of cosmic 
look-back time going from 1 to 10 Gyr. The evolution of the 
fitted scaling parameters as a function of redshift is shown 
in the left column of Fig. 1101 The corresponding evolution 
as a function of cosmic look-back time can be found in the 
right hand column. The Kormendy parameter a is shown in 
the top panels, the Faber- Jackson parameter b in the center 
panels and the FP parameters c and d in the bottom pan- 
els. Each different cosmology is represented by a different 
linestyle, listed in the insert at the top left hand frame. 



5.2 Evolution of the Faber- Jackson relation 

Evolutionary trends for the Faber-Jackson relation are com- 
parable to that seen in the Kormendy relation. No dis- 
cernible trends are found in the open cosmology, while all 
of the other high density Universes do show a mild decrease 
from b ~ 0.35 at z ~ 2 to 6 « 0.32 at z « 0. When assessing 
in terms of cosmic time (center right panel), we observe a 
near uniform increase of b over the last 8 Gyr. 

In most studied cosmologies, with the possible excep- 
tion of the ACDM02 cosmology, we find a marginal trend 
of the Fundamental Plane parameter c to decrease for z < 2, 
more or less in the past ~ 6-7 Gyr. At earlier epochs such a 
trend is entirely absent. No significant evolution of the FP 
parameter d can be observed in Fig. [TO] 



5.1 Evolution of the Kormendy relation 

For all cosmologies the evolution of the Kormendy relation 
is marginal at best. In the case of the low flm ACDM02 
cosmology we can not discern any significant change of the 
parameter a, (this may in part be due to the large uncer- 
tainties in the calculated parameter resulting from the low 
number of halos in this simulation) . In the case of the other 
cosmologies we find no noticeable change of a before a red- 
shift z ~ 2, followed by a mild increase from a ~ 0.3 to 
a « 0.38 at a ~ 0. This is also clearly visible when assessing 
the evolution in terms of cosmic time, as can be seen in the 
top right panel. 



5.3 Evolution of the Fundamental Plane 

No significant evolution has been found for the Fundamental 
Plane parameters c and d (see Fig. 1101 lower panels) . Evolu- 
tion of the Fundamental Plane mainly concerns its thickness. 

In Fig. [TT]we show the development of the FP thickness 
as a function of cosmic expansion factor aexp{t) ~ 1/(1 + z) 
for four cosmological models, and in Figs. [13] and [12] we 
show the evolution of the spread of points with the FP as a 
function of redshift in the ACDMF2 model. 

We see a systematic increase of FP thickness over the 
whole cosmic evolution in the case of the high Qm SCDM 
cosmology. While we do see a rise of the FP thickness be- 
fore a^xp < 0.5 in the ACDMF2 and ACDMC2 cosmologies, 
after that time the increase levels off and may even fiatten 
completely. Note, however, that these simulations do not 
attain sufficient halo mass resolution at higher redshifts: in 
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Figure 12. Shifting location of the cluster halo population within the Fundamental Plane. The depicted halo sample is the one in 
ACDMF2 cosmology, and is shown at four different redshifts: z = 2.61 (top left panel), z = 1.61 (top right panel), z = 0.89 (bottom left 
panel) and 2 = (bottom right panel). The abscissa and ordinate axis are arbitrarily chosen, mutually perpendicular, axes within the 
FP plane defined by (logr^, /.t, logcrt,) at 2 = fEgn. [30) l. 



these cosmologies halos still are low mass objects at these 
epochs. The one exceptional cosmology is that of the low Qm 
Universe ACDM02. Except for a rather abrupt and sudden 
jump in FP thickness at Uexp ~ 0.3, there is no noticeable 
change at later epochs. By aexp = 0.3 nearly all its clusters 
are in place and define a Fundamental Plane that does not 
undergo any further evolution. 

In summary, the trend seems to be for initial increase 
of the FP thickness followed by a convergence to a nearly 
constant value. The epoch of convergence is later for higher 
values of Qm' while the thickness remains constant for the 
low flrn ACDM02 cosmology, it involves a slow but contin- 
uous increase in the SCDM cosmology. 

On the basis of their study of galaxy merging, 
iNipoti et all l|2003l ) argued that the disposition of galax- 
ies in the Fundamental Plane is not simply a realization of 
the virial theorem, but contains additional information on 
galaxy structure and dynamics. This should be reflected in 



the location of the halo population within the Fundamental 
Plane. 

Figs. [12] and [13] show how the location of the clusters 
within the plane shifts as time proceeds. The color scheme is 
the same as for Fig. [2] Fig. ll2lshows the location of the clus- 
ters in the ACDMF2 cosmology in the Fundamental Plane 
inferred for the current epoch, ie. at redshift z = 0, 

logrh = 0.41 + 1.88 logcr„ + C/p,L . (30) 

To locate their position within the Fundamental Plane, we 
use the (artificial) coordinates Fi and F2 of the halo points 
with respect to two mutually perpendicular normalized vec- 
tors in the Fundamental Plane at z = 0, wrt. the coordinate 
system defined by the FP quantities (log r^, /i, log CTi,) (note 
that Fi and F2 do not have a specific physical significance) . 
From the panels in the figure we see that the evolution of ha- 
los involves a gradual shift along an almost universal Funda- 
mental Plane. It also shows that the halo population seems 
to evolve from a more scattered and somewhat looser one 
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Figure 13. Shifting location of tiic cluster halo population within the Fundamental Plane. The depicted halo sample is the one in SCDM 
cosmology, and is shown at four different redshifts: z = 2.61 (top left panel), z = 1.61 (top right panel), z = 0.89 (bottom left panel) 
and z = (bottom right panel). The abscissa and ordinate axis are arbitrarily chosen, mutually perpendicular, axes within the FP plane 
defined by (log r^, ^t, log cr^) at z = (Eqn.lsTJl. 



into a tightly elongated point cloud at the current epoch, 
providing interesting clues towards understanding the clus- 
ter virialization process. 

In the same vein, Fig. [13] follows the changing location 
of clusters in the SCDM cosmology in the corresponding 
Fundamental Plane at z = 0, 



log r-h = 0.37^ + 1.78 loga^, + C/p,s ■ 



(31) 



Similar to the ACDMF2 cosmology, we find that the clus- 
ter point cloud appears to assume a clearer mass stratifi- 
cation as time proceeds. While the population of clusters 
in the SCDM cosmology also appears to shift its location 
along the Fundamental Plane as it evolves, we do not find a 
trend towards a more tightly point cloud that we see in the 
LCDMF2 cosmology. We will investigate these evolutionary 
trends in more detail in an upcoming study, we have found 
indications for a possible infiuence of the different cluster 
halo merging histories in SCDM and ACDMF2 in explain- 



ing the different behaviour of the cluster point clouds in the 
Fundamental Plane. 



6 MERGING AND ACCRETION 
DEPENDENCE 

Figures Fig. [TD] and Fig. [TT] show that the evolution of the 
parameters defining these relationships is very erratic. This 
testifies to the fact that in hierarchical structure formation 
scenarios the formation and evolution of halos is hardly a 
quiescent and steadily progressing affair. Rather, halos grow 
in mass by steady accretion of matter from its surrounding 
as well through the merging with massive peers. Even the 
accretion is not a continuous and spherically symmetric pro- 
cess: most matter flows in in a strongly anisotropic fashion 
through filamentary extensions into the neighboring large 
scale matter distribution. As a result, we can expect that 
many halos will not have settled into a perfect virial state. 
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Figure 14. Accreting vs. Merging Halo Evolution; the evolution of two different halos in the ACDMF2 cosmology. Each panel concerns 
a 5h~^ Mpc comoving size box centered on the core of the halo. The sequence runs from z = 2.61 (top left panel), via z = 1.61 (top 
right panel), to z = 0.89 (bottom left panel) and finally the present epoch z = (bottom right panel). The circles indicate the location 
of HOP identified halos, with the size of the circle being proportional to the (virial) radius of the halo (overlapping circles are due to 
the projection of the corresponding spheres). Top four panels: a quiescently evolving accreting halo. Bottom four panels: a strongly 
hierarchically evolving merging halo. 
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This will certainly be the case for halos that recently suffered 
a major merger with one or more neighboring clumps. 

The detailed accretion and merging history is a func- 
tion of the underlying cosmology. Low density cosmologies 
or cosmologies with a high cosmological constant will have 
frozen their structure formation at early epochs. The halos 
that had formed by the time of that transition will have had 
ample time to settle into a perfect virialized object. Also, 
there is a dependence on the power spectrum of the corre- 
sponding structure formation scenario. Power spectra with 
a slope n < —1.5 (at cluster scales) will imply a more ho- 
mologous collapse of the cluster sized clumps, less marked 
by an incessant bombardment by smaller clumps. It may be 
clear that a more violent life history of a halo will usually 
be reflected in a substantial deviation from a perfect virial 
state. 

In order to investigate the implications of a difference 
in accretion or merging history of halos, we have split the 
samples of cluster halos in each of our cosmologies into a 
merging sample and a accretion sample. Possible differences 
in their virial state should be reflected in the quality of the 



scaling relations, in particular that of the thickness of the 
Fundamental Plane. 

The merger sample consists of those halos that suffered 
a merger with another halo that contained at least 30% of its 
mass. Fig. 1141 shows two examples of halos in the ACDMF2 
cosmology. The top sequence of 4 panels shows the evolu- 
tion of a quiesencently evolving accretion halo, by means 
of the particle distribution in a 5h~^ Mpc box (comoving 
size) around the cluster core, at z=2.61, z=1.61, z=0.89 and 
2 — 0.00. The circles indicate the location of the HOP iden- 
tified halos, with the size of the circle proportional to the 
radius of the halo (note that the overlap of circles is due to 
projection of the halo spheres). The lower group of 4 panels 
shows the particle distribution at the same redshifts for a 
halo belonging to the merging sample. Its gradual hierarchi- 
cal buildup is directly visible as the the continuous infall of 
clumps at each timestep. 

In Fig. [15] we show the evolution of the thickness of the 
Fundamental Plane for each of the two samples in the four 
indicated cosmologies. Note that our simulations do not have 
sufficient resolution for reconstructing the precise merging or 
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Figure 16. Comparison between the Fundamental Plane in the 
Lamf)daCDMF2 cosmology for clusters that underwent a major 
merger (blue dots) and clusters that followed a more quiescent 
accretion history (red dots) . The plot depicts the relation between 
harmonic radius and the quantity Y = + dloga + Cfp, in 
which c and d are the FP scaling parameters. 



accretion history before Uexp ~ 0.3 — 0.4, so that we may not 
draw conclusions on the rise of the FP thickness up to that 
epoch. Also, in the case of the ACDM02 scenario we do not 
have enough cluster halos to be able to detect any systematic 
differences between the merging and accreting halos. 

In the more recent history we do find some significant 
differences between merging and accretion-only halos in var- 
ious cosmologies, in particular the ones with a high flm. 
There does not seem to be a systematic difference between 
these groups in the ACDMF2 cosmology. The total absence 
of any difference between the Fundamental Plane of merg- 
ing and accreting cluster halos at present (Fig. I16|l . is the 
outcome of an evolutionary history without any significant 
differences between the two subsamples (Fig. 1151 lower left- 
hand panel). 

The story is quite different for the ACDMC2 and SCDM 
cosmology. While the cluster halos that undergo a major 
merger do reveal a constantly growing FP thickness, their 
accretion-only clusters do not display such a systematic in- 
crease. Instead, their FP thickness remains lower and levels 
off. In other words, accretion halos (dotted lines) do on aver- 
age display a tighter FP relation. This is particularly true at 
the current epoch. Apparently, the absence of violent mass 
gain in the case of accretion halos implies that they have 
more time to relax and virialize. This, in turn, is reflected in 
a thickness of the Fundamental Plane which does not evolve 
any further. Interestingly, it is also reflected in the radii of 
the halos (Fie:. ll7p : while the harmonic radius and half-mass 
radius of accreting cluster halos are mostly in accordance 
with each other, though with a larger spread than in the 
case of ACDMF2 clusters (lower panel), the SCDM halos 
that underwent major mergers do appear to be responsible 
for the substantial differences between the harmonic and 
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Figure 17. Comparison between the mean harmonic radius 
and the half-mass radius of the halos in the SCDM (top) and 
ACDMF2 (bottom) scenarios. The cluster samples are split into 
the clusters that underwent a major merger (blue dots) and the 
ones that accreted matter in a more quiescent fashion (red dots). 
See text for further explanation. 

half-mass radii that we see in Fig. [2l From this we conclude 
that the accretion history is a major factor in determining 
the character of the Fundamental Plane, via the impact of 
mergers on the mass distribution within halos and hence 
their radii. We discuss this in more detail in the next sec- 
tion [3 

The implications of this finding might be far-reaching. 
Given the remarkable robustness and stability of the Fun- 
damental Plane, any deviation of individual clusters from 
the FP may be a direct reflection of its recent dynamical 
evolution. This would be true if the thickness of the plane 
would be entirely due to the merger history of the clusters. 
It is certainly a viable implication of our conclusion that 
the Fundamental Plane's definition - the average plane of 
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a large sample of clusters - is nearly unassailable while we 
find strong fluctuations and deviations from the average FP 
in small samples of actively evolving clusters. 

In practice, it might mean that one could take samples 
of clusters in different redshift bands and reliably average 
them in each band to use the resulting Fundamental Plane 
to study redshift evolution of observed samples. It would 
also mean that within each redshift band you know which 
ones have had active lives. 



7 RECONCILING THE SIMULATIONS WITH 
THE VIRIAL THEOREM 

7.1 The Virial Theorem 

The Fundamental Plane is a direct reflection of the virial 
theorem which, under particular assumptions, relates the 
averaged velocity dispersion and radius of a system directly 
to its mass. All "virialised" objects will lie on a plane de- 
fined in the space of those three variables. There is not even 
any freedom in the parameters for that plane: its slope and 
location are fixed for all virialised objects. 

There are complications when assigning data to a Fun- 
damental Plane. Firstly, in its simplest form, the virial theo- 
rem assumes that the virialised objects are isolated spherical 
systems and, importantly, that they are stationary. The sys- 
tems we study are not spherical and they are certainly not 
stationary: they are generally in a state of dynamical evolu- 
tion. The possible exception to this might be the largest 
most isolated systems. Secondly, observed data does not 
have direct knowledge of the system mass except through 
interpreting the light that is observed. The universality of 
the Fundamental Plane allows us to turn the problem around 
and determine the dependence of light on mass in order that 
systems should fit on the Fundamental Plane. The simplest 
approach to this is to assume that the mass to light ratio in 
the observed waveband is directly related to mass. 

There are further issues. For example, what do we mean 
when we refer to "averages" of quantities? Using a differ- 
ent averaging process yields a different Fundamental Plane. 
There is also the fact that astrophysical systems are observed 
only in projection. 

Having said that, we can express the Virial Theorem 
in terms of the variable we have used here to describe the 
Fundamental Plane. With the notation that a virialised sys- 
tem of mass M has a velocity dispersion V , half mass radius 
Thaif and harmonic radius we have, up to normalising 
constants: 



M 



M 



' half 



(32) 



where E is the projected (surface) mass density. Eliminating 
M from these and taking logs yields an expression for the 
Fundamental Plane: 



logrh = 



Th 
Thalf 



^ 2\ogV + QA^ji 



(33) 



where we have transformed the surface mass density S into 
logarithmic astro-units via 

/I = -2.5 log E (34) 
We have explicitly written equation (|33p in such a way as 



to expose the different roles of harmonic and half-mass (ge- 
ometric) radii. The relationship between these radii in our 
models is illustrated in Fig.[2]and in Fig. 1171 the latter differ- 
entiating between merging and quiescently accreting halos. 



7.2 Renormalising the FP Simulations 

It is important to understand why the coefficients of the 
model Fundamental Plane might differ from the expecta- 
tions based on the use of the virial theorem. Luminosity is 
not involved here so we cannot appeal to a varying mass-to- 
light ratio. Moreover, the model Fundamental Plane is well 
defined and so we cannot say that this is merely a question 
of fitting. 

The are at least two possible sources for this systematic 
difference between the model and the virial theorem. The 
first is to blame the HOP technique and assert that it sys- 
tematically underestimates the cluster masses. The second 
is to say that the internal cluster properties (like velocity 
distribution) vary systematically with mass an so the nor- 
malisation of the virial plane is mass dependent. 

Either way, we shall model in a mass dependency and 
consider this in relation to the HOP technique. The process 
for the variable virial normalisation is analogous. 

The samples of clusters derived from these simulations 
are all based on the HOP technique. There may well be 
a systematic bias in the assignment of particles to clusters 
(see section [5TT| . As a consequence, the radii and velocity 
dispersion derived for a HOP selected cluster will also be 
biased. Clearly the bias will be more significant for smaller 
systems. 

In this subsection we seek to account for systemic ef- 
fects of using HOP for identifying cluster membership, and 
derive a renormalisation procedure taking account of this 
and matching the dataset to the expected virial theorem 
Fundamental Plane (equation [33} . 

The easiest way to model this bias is to assume that the 
model-based estimate (biased) for the mass, M, is related 
to the actual mass by a simple scaling relationship 

M 

— cxA^— (35) 

for some exponent a. The virial expression for the mass then 
becomes 
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where the right hand side now refers to quantities derived 
from the model. We can eliminate M from this in terms of 
the model surface mass density E = M/rf^ to give 
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Taking logs and using ^ 
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(38) 



which is the expression for the Fundamental Plane in terms 
of the (biased) model derived quantities. This should be 
compared with equation (|33p : we see how the bias modelled 
by a affects the position and slope of the virial Fundamental 
Plane. 
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Model 


d 


a 


c 


c(modcl) 


SCDM 


1.78 ±0.01 


0.058 ± 0.003 


0.38 


0.37 ±0.031 


LCMDFl 


1.69 ±0.084 


0.084 ± 0.024 


0.37 


0.38 ±0.016 


LCDMF2 


1.88 ±0.04 


0.031 ± 0.011 


0.39 


0.41 ±0.011 


LCDMF3 


1.81 ±0.02 


0.050 ±0.014 


0.38 


0.38 ±0.056 


OCDMOl 


1.60 ±0.083 


0.111 ± 0.024 


0.36 


0.35 ±0.017 



Table 5. Fitting biased models to ideal virial Fundamental Plane. 
The parameter a emulates the limitations of the HOP group 
finder. 



The procedure now, for each simulation, is to select a 
value of a that makes the coefficient of log V in equation (|38p 
equal to the virial value 2. That a then allows a calculation 
of the coefficient of /i that can be compared with the value 
derived from the simulation. The results for a selection of 
models is shown in Table [5] 

The conclusion to be drawn from this is that, for each 
model, there is indeed a value of the a parameter that re- 
produces the Fundamental Plane fits for the models. 

T.3 Observed cluster FP 

The best availabl e data set is the ENACS data of 
l|Adami et al.lll99^ . Equation IfTT} describing that Funda- 
mental Plane, in the current notation, reduces to 



logiie = (0.49 ±0.05)^ + (1.12 ±0.11) log (j„ 



(39) 



There is considerable uncertainty in this relationship: the 
coefficient of log a is quite far from the ideal 2.0 and the 
coefficient of ^ is higher than the nominal 0.4. 

The usual way to reconcile this with the virial Funda- 
mental Plane is to argue that the mass to fight ratio of the 
cluster sample is mass dependent: 

^ oc Af" (40) 

Using an argument that parallels the derivation of equation 
(|38p . the Fundamental Plane expressed in terms of velocity 
and surface mass density is 



log rh 
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(41) 



The data give 13 — 0.28 ± 0.19 which gives rise to c = 
0.31 ± 0.02, a long way from the data-derived 0.49. It is 
clearly not possible to reconcile the ENACS data with the 
virial theorem Fundamental Plane, let alone the numerical 
simulations. 



8 CONCLUSIONS AND DISCUSSION 

We have studied three structural scaling relations of galaxy 
clusters in thirteen cosmological models. These relations are 
the Kormendy relation, the Faber- Jackson relation and the 
Fundamental Plane. Their validity and behavior in the dif- 
ferent cosmological models should provide information on 
the general virial status of the cluster halo population. The 
cosmological models that we studied involved a set of open, 
flat and closed Universes with a range of matter density pa- 
rameter Qrn and cosmological constant flA- 



The cluster samples are obtained from a set of A'^-body 
simulations in each of the cosmologies. These simulations 
concerns a box of 200/i~^Mpc with 256^ dark matter parti- 
cles. The initial conditions were set up such that the phases 
of the Fourier components of the primordial density field are 
the same for all simulations. In this way, we have simulations 
of a comparable morphological character: the same objects 
can be recognized in each of the different simulations (be it 
at a different stage of development). 

After running the simulations from z = Ato the current 
epoch using the GADGET2 code, we used HOP to identify 
the cluster halos. We investigated whether each halo popu- 
lation obeyed a mass-radius relation akin to the Kormendy 
relation, a mass-velocity dispersion relation similar to the 
Faber- Jackson relation and a two parameter family between 
mass, radius and velocity dispersion that resembles a Fun- 
damental Plane relation. We studied the dependence of the 
obtained scaling parameters as a function of the underlying 
cosmology and investigated their evolution in time. 
Our results can be summarized as follows: 

• In each cosmological model we do recover Kormendy, 
Faber- Jackson and Fundamental Plane relations for the pop- 
ulation of cluster halos. This is a strong indication that the 
halos are in a virialized state, as expected in hierarchical 
clustering scenarios. 

• There are significant differences between the measured 
parameters of the various scaling relations and those seen in 
the observational data. Our fit for the FP in the ACMDF2 
model is 



log Th = 0.41/i + 1.86 log (7 + const. 



(42) 



This can be reconciled with the expectation from the virial 
theorem, but not with the ENACS Fundamental Plane. 

• We do not find any significant dependence of the pa- 
rameters a and h of the Kormendy and Faber- Jackson rela- 
tions on the value of fim- There is also no indication for any 
influence of Q,a on the scahng relations. 

• While the FP parameters c and d axe not dependent on 
Q.m and Q,a, there is a slight suggestion that the Fundamen- 
tal Plane would have a lower thickness for low Q,m ~ 0.1 
cosmologies. 

• With the exception of low Q.m Universes, we find a mild 
increase of the Kormendy parameter a and a mild decrease 
of the Faber-Jackson parameter b from z = 1 to the present 
epoch. From 2 = 4 to z = 1 we did not find any discernable 
evolution. 

• While the Fundamental Plane parameters c and d do 
in general not show a signiflcant evolution, the higher Q.m 
cosmologies do involve a slight decrease of FP parameter d 
during most recent epochs {z <2). 

• The thickness of the Fundamental Plane does evolve sig- 
nificantly, with an initial increase followed by a convergence 
to a more ore less constant value. The convergence epoch is 
later for higher density cosmologies. This probably reflects 
the gradually virializing tendency of the cluster population. 

• Given our expectation that there is a difference in virial 
state between quiescently accreting clusters and those ex- 
periencing massive mergers, we have investigated the evo- 
lution of the Fundamental Plane thickness for samples of 
merging clusters and samples of accreting clusters. We find 
that accreting clusters at recent epochs do appear to be bet- 
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ter virializod than the merging population and that the FP 
thickness is smaller in the former. 

• We find that for all investigated cosmologies the Fun- 
damental Plane is remarkably stable, despite the enormous 
evolution of the individual systems. The only significant evo- 
lution, that of its thickness, might be due in a large part to 
the importance of merging of individual systems. 

• If indeed the thickness of the Fundamental Plane might 
be entirely due to the merger history of the cluster halos, the 
distance of an individual cluster to the Fundamental Plane 
would be a direct refiection of the cluster history. 

• We see direct evidence that major mergers have effected 
the relationship between the galaxy haloes in the cluster in 
that the relationship between the half-mass and harmonic 
radii is disturbed. Nonetheless, the evidence from the models 
tells us that this does not affect the slope of the Fundamental 
Plane: clusters that have undergone major mergers lie in the 
same place as those that have grown by steady accretion. 

Finally, what is desperately needed is better data on 
the cluster fundamental plane. We might speculate that the 
distance of a cluster from the plane defined by the data 
somehow reflects the evolution of the cluster, but we will 
not get evidence for the hypotheses derived from numerical 
experiment until there is more high quality data. 
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